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Abstract 

One-loop corrections to the helicity amplitudes of all 2 ^ 2 subprocesses are calcu- 
lated in QCD and in N=l supersymmetric Yang- Mills theory using two versions of 
dimensional regularization: the 't Hooft-Veltman scheme and dimensional reduction. 
Studying the structure of the soft and collinear singularities, we found universal tran- 
sition rules for the squared matrix element which can be used to translate the results 
obtained in these schemes to the results valid in the conventional dimensional regular- 
ization scheme. With explicit calculation it is demonstrated that the one-loop helicity 
amplitudes of the 2 — > 2 subprocesses calculated using dimensional reduction in the 
N=l supersymmetric SU{N) gauge theory respect the supersymmetry Ward identities. 
Our transition rules can also be used to calculate the next-to-leading order Altarelli- 
Parisi kernels in the dimensional reduction scheme when they satisfy supersymmetry 
Ward identities as well. 



Work supported in part by the Schweizerischer Nationalfonds 



1 Introduction 



This paper discusses the consistency of various versions of the dimensional regularization 
schemes for regulating both ultraviolet and infrared divergences in calculations of infrared 
safe quantities for hadron-hadron collisions at next-to-leading order in QCD and in N=l 
super symmetric Yang-Mills theory. The discussion is based upon the calculation of the one- 
loop radiative corrections to the helicity amplitudes of all 2 2 parton scattering processes. 
The motivation for this study is a recent paper by Bern and Kosower in which the authors 
worked out string-theory-based technique for evaluating multi-gluon amplitudes in one-loop 
order. The power of the method is demonstrated clearly by the recent calculation of the one- 
loop corrections to five gluon helicity amplitudes using the string-theory-based technology 
[0]. The field theory interpretation of the new results suggested that the use of background 
Feynman gauge, the helicity method and a new version of dimensional regularization, the 
so called four- dimensional helicity scheme, results in great technical advantages also in the 
standard field theory calculations. 

In perturbative QCD, with the evaluation of the loop corrections one does not obtain 
yet physical cross-sections. The Bremsstrahlung contributions have to be added as well. 
In particular, finite hard scattering cross-sections are obtained only (after trivial ultraviolet 
renormalization) by cancelling soft and coUinear singular terms between loop corrections and 
Bremsstrahlung contributions, and by subtracting the initial-state collinear singularities. In 
ref. 0, the singular terms appearing in this procedure have been worked out analytically 
using the conventional dimensional regularization scheme. The analytic expression obtained 
in ref. for the soft and collinear contributions have universal, process independent charac- 
ter, but they depend on the regularization scheme used. The sensitivity to the regularization 
scheme, however, are also exhibited analytically, therefore, the result of ref. can easily be 
transformed to the various versions of the dimensional regularization schemes. 

Physical cross sections of infrared safe quantities are obtained by folding the finite hard- 
scattering cross section with parton densities. In a next-to-leading order calculation the 
Q^-evolution of the parton densities has to be carried out with the next-to-leading order 
Altarelli-Parisi kernels which also depend on the regularization, factorization and ultraviolet- 
subtraction schemes. In the conventional MS factorization scheme the Altarelli-Parisi kernels 
have been calculated by Curci, Furmanski and Petronzio [Q. Changing the regularization 
scheme will generate changes in the next-to-leading order Altarelli-Parisi kernels, in the 
hard-scattering cross sections and in the value of Aqcd, but the physical cross section has 
to remain the same. Therefore, calculating the hard-scattering cross sections in different 
schemes, we can identify the universal transition terms which define the transformation of 
the next-to-leading order Altarelli-Parisi kernel from one scheme to another. 

Recently, it has become clear that with appropriately modified dimensional regularization 
schemes the advantages of the helicity method can be maintained also for loop calculations 
[H, § . Three modifications have been proposed: the 't Hooft- Veltman scheme, regularization 
by dimensional reduction and the so called four- dimensional helicity schemeQ. Changing 
the regularization scheme, the soft and collinear singular contributions will change but this 

^The four-dimensional helicity scheme was proposed in the context of the string-theory derivation of the 
one-loop corrections to the four-gluon helicity amplitudes According to our knowledge, no operational 
definition of this regularization has been defined in the standard field theory framework. 
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change has universal character. Therefore, it will be sufficient to find the changes generated 
in the hard scattering cross sections only for few processes. Then one knows the transition 
functions for all other processes. Deriving the transition terms has great practical importance 
because the phenomenological fit of the parton densities and their Q^-evolution is worked 
out in the conventional regularization scheme. Therefore, if we calculate the next-to-leading 
order corrections in some other schemes, we must be able to transform the result into the 
conventional scheme. 

In this paper, with the explicit calculation of the next-to-leading order corrections to all 
2-^2 parton-parton scattering amplitudes, we determine the transition rules necessary to 
translate the results obtained in the 't Hooft-Veltman and dimensional reduction schemes into 
those in the conventional scheme. The rules are process independent, therefore, they can be 
applied to other processes. For example, calculating the next-to-leading order corrections to 
the 2^3 scattering amphtudes, one can use such a regularization scheme that is consistent 
with the application of the helicity method and supersymmetry (dimensional reduction) and 
then one can easily translate the result into the phenomenologically relevant conventional 
scheme. 

It has been noted long time ago that in calculating scattering amplitudes of massless 
partons in perturbative QCD, supersymmetry can be used as a technical tool to check and 
simplify the calculation. The reason for this is that if we change the color representation 
of quarks to the adjoint representation of a Majorana fermion we obtain the N=l super- 
symmetric Yang-Mills theory in the Wess-Zumino gauge. If we use dimensional reduction 
which respects supersymmetry, then we can test supersymmetry Ward identities between the 
amplitudes of various subprocesses with different number of external fermions and bosons. 
These relations, which are valid in all order in perturbation theory, provide us a very signif- 
icant check on the correctness of the calculation. One can also turn the argument around 
and use the supersymmetry Ward identities to obtain gluonic amplitudes from the fermionic 
ones correcting for the differences in the contribution of fermion loops. 

We carried out all the calculations presented in this paper using standard Feynman 
diagram method, but using all the advantages of the helicity method, the background Feyn- 
man gauge and the 't Hooft-Veltman and dimensional reduction regularization schemes. The 
string theory derivation of the gluonic amplitudes in ref. 0] suggests that significant technical 
advantages can be achieved in addition with cleverly combining class of Feynman diagrams 
and reducing tensor integrals to scalar integrals. We did not make attempts to simplify our 
calculation in these directions. However, we used color subamplitudes and the method 
of ref. [0] to reduce the tensor loop integrals to scalar integrals which is easily adaptable to 
calculations in massless QCD. This part of our field theory calculation still cannot match 
the simplicity achieved by the string-based technique. 

The organization of our paper is as follows. In sect. 2, we describe the various versions of 
dimensional regularization briefiy. In sects. 3 and 4, we give an overview of the techniques 
which enter our calculations: the helicity method and color decomposition of the amplitudes, 
the reduction of tensor integrals and use of background field method. These technicalities 
have been discussed in much details elsewhere in the literature, therefore, our presentation 
is constrained to setting the notation. Sect. 5 contains the one-loop color subamplitudes in 
QCD which constitutes one of our main results. We present the results in the unphysical 
region where all kinematic variables are negative which may become useful in the future. We 
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also calculate the shift in the strong coupling constant caused by using different regularization 
schemes. In sect. 6, we calculate the next-to- leading order loop-correction to the square of 
the matrix element in order that we can compare our results to that of Ellis and Sexton 
[§]. We establish transition rules among the one-loop amplitudes obtained using different 
versions of dimensional regularization. Sect. 7 presents the one- loop color subamplitudes in 
the supersymmetric limit of QCD, N=l supersymmetric Yang-Mills theory, and elaborates 
on the usefulness of supersymmetry in QCD calculations. In order to find the complete 
transition rules among the various regularization schemes, in sect. 8, we discuss the difference 
among the Bremsstrahlung contributions obtained using different regulators. In sect. 9, we 
shall spell out how the hard-scattering cross section changes with changing the regularization 
scheme. In sect. 10, we shall discuss how the change in the hard-scattering cross section due 
to the use of different regularization schemes can be compensated by a change in the parton 
density functions. In this way, we can present some consistency checks on the correctness 
and process independent property of our transition rules. Finally, sect. 11 will summarize 
the main results. 



2 Regularization schemes within dimensional regular- 
ization 

In massless QCD, the conventional dimensional regularization is an attractive scheme because 
it simultaneously regulates infrared and ultraviolet divergences, is manifestly gauge invariant, 
consistent with unitarity and simple to implement. 

When defining the different versions of dimensional regularization both for loop and phase 
space integrals, it is useful to distinguish two class of particles: observed and unobserved 
ones. Unobserved particles are those virtual ones which circulate in internal loops as well as 
those which are external but soft or coUinear with other external particles. All the rest are 
observed particles. Unitarity demands that unobserved particles are treated uniformly. 

The most important ingredient of dimensional regularization is the continuation of the 
momenta of the unobserved particles into d ^ A dimensions, thus rendering the integrals 
over these momenta finite. Having this done, one is left with a lot of freedom how to treat 
the momenta of observed particles and the polarization vectors of all particles. The different 
choices lead to different versions of dimensional regularization. 

The original choice made by 't Hooft-Veltman [y] was to continue both the momenta 
and the helicities of unobserved particles into d dimensions, while leaving the momenta and 
helicities of observed particles in four dimensions.^ 



There is a third version of dimensional regularization, the dimensional reduction |Tl 
widely used in connection with supersymmetric theories. This technique consist of continuing 
the virtual momenta of loop integrals into dimensions < 4, but keeping all polarization 
vectors in four dimensions. Because of d < 4, the method is not directly applicable for 
regulating infrared divergent phase space integrals. It is essential that < 4 because the 



^As a matter of fact, the original proposal was meant for ultraviolet regularization of loop integrals only. 
The possibility of using it also for regulating infrared singularities was initiated in ref. [nO|. For a more 



detailed recent application see ref. ||], 
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four-dimensional vector field is split into a (i-dimensional vector plus a field which transforms 



as (4— (i)-dimensional scalars under gauge transformation [12|. It can be shown, however, that 



the use of this Lagrangian is equivalent to using the original Lagrangian without the splitting 
of the gauge field and perform the algebra in four dimensions. The only subtle point which is 
a reminiscent of the splitting is that one has to distinguish between four-dimensional metric 
tensors coming from the Feynman rules and d-dimensional ones emerging from momentum 
integrals with more than one loop momentum in the numerator. In this way, gauge invariance 
is maintained which has been checked explicitly up to two loops [1^. We remark that the 
operational definition of dimensional reduction is not complete in complicated cases (see e.g., 
refs. [|1^, Q). However, this incompleteness does not affect our calculation. 

In supersymmetric theories, there is a different formulation of dimensional reduction 
when supergraph technique is applied. In this approach all algebra is performed in four 
dimensions until only final scalar integrals are left in which the continuation of the loop 
momenta is performed. In ref. [Q, for relatively simple calculations the equivalence of the 
two approaches has been pointed out. 

Inspired by string-based rules for calculating loop amplitudes in gauge theories, in ref. 

a new version of dimensional regularization, the four-dimensional helicity scheme has 
been proposed. The operational definition in the string calculation is such that only the 
momenta of unobserved particles are continued into d dimensions, all polarization vectors 
and the momenta of observed particles are kept in four dimensions. The meaning of such a 
definition in field theory is obvious in most of the algebra except for the reduction of tensor 
loop integrals. In the string-based calculation, the loop momentum is integrated out at the 
string level which implies that the reduction of tensor integrals has been implicitly achieved 
in four dimensions. Therefore, in the field theory version of the four- dimensional helicity 
scheme, one may try to set li = 4 in the algebraic part of tensor loop integrals and treat the 
remaining scalar loop integral in d dimensions. In the present paper, we do not pursue this 
possibility any further. 

The requirements that a certain regularization scheme is a consistent regulator in a gauge 
field theory are that it has to respect gauge invariance and unitarity. Gauge invariance is 
known to be maintained by the conventional, 't Hooft-Veltman and dimensional reduction 
schemes. Unitarity is obviously maintained if the momenta and helicities of all particles are 
continued into d dimensions, therefore, the conventional scheme is considered a consistent 
scheme in massless QCD. The other three versions of dimensional regularization do not treat 
momenta and helicities equally, therefore, it is not obvious that unitarity is respected by 
these schemes. 

In axial gauge, the collinear singularities come from self-energy contributions on external 
lines only [Q, therefore, it is expected that unitarity is preserved in the 't Hooft-Veltman 
scheme. In the case of dimensional reduction and four-dimensional helicity scheme, we do 
not know such a general argument. 



3 The helicity method and color partial amplitudes 



We use the helicity method in its simplest ||T5[, crossing symmetric version |[T^. In order to 
fix our notation, we present the definitions for the spinor algebra here. For more details, the 
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reader is referred to Appendix A of ref. |T8| . 

Let V'(p) be a four-dimensional spinor satisfying the massless Dirac equation: 

=0,p^ = 0. (3.1) 

We define the two hehcity states of ip{p) by the two chiral projections 

= 1(1 ± 75)^(p) = ii^^ip)r- (3.2) 

The second equation is a conventional choice of relative phase between opposite helicity 
spinors fixed by the properties under charge conjugation C, 

{^{p)r = C{^{p)r, C%C-' = 7^. (3.3) 



Following ref. [0, 0, we introduce a new notation 



b±) =^±(p), (p±| =^±(p), (3.4) 

(pq) = {p- \p+), [pq] = {p+ \p-). (3.5) 
With this notation the normalization of the spinor is expressed as 

{p±\-f,\p±) = 2p„ (3.6) 

and we have the useful property 

|p±)(p± 1 = ^(1 ±75)/ (3.7) 
The polarization vector of an outgoing massless vector of momentum p is defined as 

eHp,k) = ±!^^^hm (3.8) 

where k is an arbitrary (reference) momentum. Then we find that the usual requirements 
for the polarization vectors, 

e+ip,k)-e+ip,k) = 0, (3.9) 
e-^{p,k) ■e'{p,k) = -1 (3.10) 

are satisfied. 

As mentioned in the introduction, we also make use of the color subamplitudes P] which 
give gauge invariant decomposition in color space. For the tree and one-loop four-point 
amplitudes we findQ: 



■^We use the same notation as in ref. j|] to distinguish the different processes and take all particles 
outgoing. We note that in the case of gluon-gluon scattering, our decomposition differs formally from that 
given in ref. fit]. The difference is the omission of the term proportional to TrT" = 0. 
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AT{q,Q;Q,q) 
AT'iq, q; q.q) 

AT\g.9:qA) 



g ( '^ni3'^i2i4 ~ ^'^iii4^j2i3 ) '^4;o(l) 2; 3, 4) 



1 



9^ X! \^iM^)^i2o(A) - -^^na(4)'^i2<7{3) 

x64;o(l,2;a(3),(T(4)). 



g2 j2 (T'^^WTM.)) c4;o(^(l),a(2);3,4). 



(TGS2 



3<t{4) ^ 



(TGS4/Z4 

xci4;o (a(l),(T(2),(T(3),a(4)) 



^n«3^«2«4 jY ^iii4'^«2i3j '^4;i(l; 2; 3, 4) 

+5ili3'^i2i4«4;2(l,2;3,4)] . 



Al ^°°^{g,g;q,q) 



x64;i(l,2;a(l),(T(2)) 

+5na(3)5i2<T(4)&4;2(l, 2; (t(3), 0-(4)) 



^ (TMi)TM.))^^^^C4;i(a(l),a(2);3,4) 

f7g52 

+'^i3i4'^aia2C4;2(l,2;3,4)] . 



(3.11) 

(3.12) 

(3.13) 
(3.14) 



(3.15) 



•AI ^°°^{q,q;q,q) = g^ (SiM3)Si2a{A) - ^S,^aii)Si^a{3)) (3.16) 



(3.17) 



Al~^"°^ig,g,g,g) = g^[ ^ Tr(T''-(i)T"-(2)T'^.(3)2-a.(4)) (^318) 

o-G54/.^4 

Xrf4;l(cr(l),a(2),Or(3),Or(4)) 

+ ^ Tr(T'''"(i)T"'"(2))'jr(T'^''(=')T'*'"W) 
(TeS4/z| 

X rf4;2(a(l),(T(2);a(3),(T(4))]. 

In these equations, Sn is the permutation group of n elements, is the cychc group of n 
elements, i.e., Sn/Zn means the group of permutations of n elements with cyclic permutations 
removed. 

The hehcity subamplitudes for ra-gluon processes enjoy a number of useful properties |jl| 
which simphfy the actual calculation (or alternatively, can be used as checks of the results). 
The amplitudes of the other processes posses much less symmetry. However, useful relations 
can be discovered in the supersymmetric limit (see sect. 7). 

The Feynman diagrams for quark-quark, quark-gluon and gluon-gluon scattering have 
been presented in ref. |H]; we do not repeat them in our paper. The diagrams that enter 
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the calculation of a given partial amplitude can be found by imposing the condition that 
they have the proper color structure. In the sect. 4, we describe how to compute a one-loop 
Feynman diagram. The external legs of the diagrams to be computed are on shell, i.e., 
for massless partons the momentum squared of the external legs is zero. In dimensional 
regularization a loop insertion on a particle line is proportional to the momentum squared 
of the line. This means that those diagrams which have a loop insertion on an external line 
are zero in dimensional regularization. Motivated by the results of ref. [19|, we shall use 
background Feynman gauge. The necessary Feynman rules can be found in the standard 
reference [pH . 



4 Calculation of one-loop Feynman diagrams 

One of the methods we used for calculating one-particle irreducible Feynman diagrams is as 
described in ref. In this section, we recall the basic ingredients of the method and give 
all necessary formulae for calculating four-parton processes explicitly. 

The general algebraic expression for a one-loop one-particle irreducible Feynman diagram 
with N propagators in d = 2uj {uj = 2 — e) dimensions is 

^(M) (41) 

y (2,1)2" nf.i(9?- m? + i>))' C'l; 

where i is the loop momentum, pi are the external outgoing momenta, g,. are the momenta of 
the propagators and N{q,p) is the numerator which receives its contributions from vertices 
and numerators of propagators and contains the dimensionless coupling.^ After Feynman 
parametrization the loop integral can be performed to give 

Cutting two of the propagators in the loop produces two tree diagrams — a two-tree. A 
double cut can be performed (^^^ different ways. The set of two-trees will be denoted by 
T2. With this notation and for massless particles in the propagators 



D{x,p) = -J2iVc,) [Y.Pc] 

T2 \ c / 



(4.3) 



where Xc,, Xc^ are the Feynman parameters of the two cut lines and the sum over c is the 
sum of the external momenta belonging to one of the trees. In the numerator in eq. ( [4.2|) , 
Nq{J,p) = N{J,p) and 




Jr= E^cJEPc . (4.4) 



''^A factor of /i^^^^^^ is omitted. When a cross section is calculated this factor appears in the corresponding 
tree diagrams, as well. Therefore, for the cross section it means an irrelevant overall factor which disappears 
at the end of the calculation when physical dimension is considered. 
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In this equation one of the cut hues is r, while is the Feynman parameter of the other 
cut hne. In the sum over c, the sum of the external momenta belonging to the tree to which 
Qr flows is to be taken. Nk{J,p) for A; > is obtained from Nq{J,p) by contracting k pairs of 
J's in all possible waysQ and summing over all such contractions. If no contraction is possible 
then = 0. Contraction of J's means 



JrJg 



-,9(2^)- 



(4.5) 



This rule is valid in this simple form only if the propagators in the loop constitute a continuos 
flow of the momentum which for a one-loop diagram can always be arranged. 

We give formula (4^) for the cases that enter our calculation (A^ = 2, 3, 4) together with 
the necessary integrals explicitly in Appendix A. 

As a second method, we calculated the one-loop tensor integrals using the standard 
Passarino-Veltman type scheme . 

We note that the tree level diagrams with one-loop renormalization insertions have also to 
be taken into account. The MS ultraviolet counterterm for the n-point scattering amplitude 
at the one-loop level is well known: 



(n-2) 



-Hi 



-A 



tree / 



(4.6) 



where /3o = (llA^c — 2A^/)/3 and A^"^^ is the four-dimensional tree-level amplitude. 



5 One-loop color subamplitudes in QCD 

In this section, we present the results obtained for the QCD color subamplitudes. First, we 
list the amplitudes in the 't Hooft-Veltman scheme. Next, we give the differences of the 
amplitudes obtained in the 't Hooft-Veltman and dimensional reduction schemes. Finally, 
we calculate for dimensional reduction in terms of Aj^ of conventional dimensional 
regularization.0 

In order to check the following results, we calculated the nontrivial integrals in two 
different ways (one is described in Appendix A) and checked gauge invariance in all external 
gluon legs using longitudinal gluons. The amplitudes for gluon-gluon scattering without the 
fermionic contribution have already been calculated in ref. Our results are in complete 
agreement with the published ones. 

5.1 Results in the 't Hooft-Veltman scheme 

In the following, we constrain the presentation of the results to the minimal information 
from which after trivial algebra — i.e., using parity transformation or the cyclic property of 
four-gluon amphtudes — any helicity amplitude can be obtained. To calculate the square of 
the matrix element for process B, the same partial amplitudes are to be used as for process 

^This contraction is done in 2uj dimensions both in 't Hooft-Veltman and dimensional reduction schemes. 
^Attf- in the 't Hooft-Veltman and conventional schemes are identical. 
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A, only the color and helicity sums differ. Therefore, the subamplitudes for processes A and 
B are not written separately. The d4^-2 subamplitudes will not play any role in the rest of this 
paper. They can be obtained from the from the Nj independent part of the d^-i amplitudes 
using the decoupling equations of ref. [|l|, where they were published, therefore, we do not 
repeat them here. We remind the reader that according to our convention all particles are 
outgoing.[] 

We start our description of the results with those amplitudes which are finite in four 
dimensions (the corresponding tree amplitudes are vanishing). 

Process A: 



The amplitudes which vanish at tree level vanish at one-loop level as well. 
Process C: 



-4:1 I 5 



-;+,-) 



i (34) [12] [13] [23] 
487r2 S12S14 

i (12) (14) (24) [34] 



487r2 



2 Nr 



3N^ 



2 Nr 



S12 

S12 



<2 (+,+;+,-) 



-'4;2 ^ 



-,-;+,-) =0. 



(5.1) 

(5.2) 
(5.3) 



Process D: 



= 3^(*-^'/)(l5)(|||3)(JI)- 

i f24P 
4,1V , , , ; ^g^2 V -^^12] (23) (34) [41]^ ' ^ ' 

(5.6) 

All the other partial amplitudes are divergent in four dimensions. We present the results 
in the unphysical region, where all kinematic variables are negative, therefore, the arguments 
of the logarithms are positive. The analytic continuation to the physical region can easily 
be performed keeping in mind that the branch cut of the logarithm is obtained by inserting 
the it] associated with each kinematic variable: 

(-%)'' ^ |s,,|e--^e(^-), (5.7) 
log(-Sij) -> \og\sij\+inQ{sij), (5.8) 

where 9(x) is the usual step function. The definition of Pq, appearing in the following results, 
is given after eq. (4.6). One finds: 

^We have checked the results for Processes A and C in the 't Hooft-Vehman scheme using Passarino- 
Veltman type reduction of tensor integrals. 
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Process A: 



af^i-, -; +, +) = cra4;o(-, -; +, +)F^;i {s, S12, S13, su), 



Fa;l •^IS? ^13, S14) — 

' '2 3 11 2 



-^1 iN. 



_ , Si4 13 2 

---- + 7^ log— + — +7r2 

e 3e e S12 9 



2 10 
'3e ~ Y 



1 



2 3 2 S12 lSi4 / Si3\ / 2 "5l4 , 2 

— log 8 + 1 log + TT 

e'' e e S13 2 S12 V S12/ V S12 

+ — log — \ - -(5o. 



+; +) = cra4;o(-, +; -, +)-F S12, S13, su), 



-^a;l^(^>Sl2,Sl3, S14) = 



Sl4/ 



1 



2 3 11 2, si4 13 2 
■J5-J+3j-jl<'g- + y + -^ 

"2 3 2, S12 
— log 

. 6'^ e e Si3 



2 lo- 
ss ~ ¥. 



AT. 



1 -514 A _ S12 

2 Si3 V Si3 



log h TT H lo| 



Sl3 



Sl3 512 



1 



«?2 (-> -; +) = cra4;o(-, -; +, +)-^a;2 S12, Sl3, su), 



-^a;2 "512, Sl3, Su) 

V S14/ A^c 



2 £12 _|_ -l Si4 
£ Si3 2 S12 



--log- + -- ( 1 - - ) ( log^ ^ +7r^) + ^log^ 



Sl3 
•§12 



Sl3 



S12 Si3 



a4;2 (-> +; +) = cra4;o(-, +; -, +)F^/(£, S12, S13, su), 



-^a;2''"(^) *12, •5l3, ^14) 

I S14/ AT, 



2, S12 IS14A -512 \ A 2-514 , 2\ -514 , Si4 
— log 1 log + TT log 

e si3 2si3 V S13/ V S12 / si3 S12. 
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Process C: 



+; +, -) = crC4;o(-, +; +, -)F^./{e, S12, S13, S14), 



^c;l~i^i '^12? -513? -514) — 

■2^^ ^ ^3 3 2 



-^1 iN. 
S12 , 



- , - 1 --Lo ■ , - 2 1 1 2 -513 ,3 Si3 
^ - + - log 7^ + 7^71" - 7^ log h - log 



Sl3 _ 7 1 

£2 2e ' £ '^"^ S12 2 2' 

1^13 
2si4 



, -Jg — + — M 

Sl4 S12/ S12 VSi4/ 



S12 2 S12 
2 



2s 



1^12 
2si4 



1 log + log + TT 

V Sl4 S12/ S12 \Su/ 



1 

-/ 

e 



C4;l ( + , +, -) = CrC4;o( + , -; +, -)F^1 {£, S12, Sl3, S14), 



-^ctl "^12, -513, S14) 



Sl3 



1 S12 A 2 -513 



— - — + -log— -3 + 7r^ + - — log^ — + TT 



2£ £ S12 



2 Si4 



i^flog^^ ■ 
2si4 V S12 



S12 

-f3o. 
e 



(14) 



<2 +; +, -) = cr(-i)|^[23][24]F-+(£, s^, Sis, S14), 



-^c;2'^(^> •S12, Sl3; ^14) 



S12 \ Si2 Si3 Si3 V 513/ S12 



,2 \ s 



+ 



2si4 



1 2 •Sis , 2 
log hTT 

S12 



c^I (+, -; +, -) = cri||[13][14]F+-(£, 5i2, 513, sm), 



-^c!2 -^12) -513) -Su) 



S12 V "^12 / "5l3 "S14 V Si4 / Si3 



,2 \ s 



2s 



13 



log^ ^ + TT^' 
S12 
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Process D: 



~) +) +) — Cr'i4;o(~) ~) +) -5125 Sis, ^u) 



11 2, si4 67 
log 

3£: e Si2 9 



TT 



2 10 



(5.25) 

(5.26) 



e 



diji-, +, -, +) = Cr4;0(-, +, -, +)F^;i"^(e, Si2, Si3, s 



14 



(5.27) 



^d;li^^ -^IS; Sia, S14) — 

■2^^ ^ ^4 11 2, si4 67 

72 



(5.2^ 



Si4 



log \- 71' 

Se e S12 9 



'13 



Si4 



2 10 

3£ ^ Y 



1 - ( ^ - ^ log 

■S13 -Sis/ -512 

— ^ - 2 (log — + ^ 

, -Sis / \ ■^12 



+ 



11 



Nr - -rN^ 



Si4 , Si4 3 S12S 

— log Nf 

Sl3 ^Si2 2 ^ 
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, ,log^^+7r^U--/5o. 

'13 V S12 



In these equations 



1 r^\-e)r{l + e) 



is a ubiquitous prefactor, = A*"^ — 1. For the sake of completeness, we list the tree level 
expressions here: 



o-i-fii—, —', +, +) 

a4;o(-,+; -,+) 

c4;o(^,^;g^,g~) 



. (12) [34] 

Sl4 

.(13) [24] 



Sl4 



{Iq){Iq) 



^\3 



(12)(2g)(gg)(gl)' 
■ (12) (23) (34) (41)- 



(5.30) 
(5.31) 

(5.32) 
(5.33) 



For process C, c^-q is vanishing if the gluons have the same helicity. In eqs. ( ^.32|) I denotes 
the negative helicity gluon. For process D, d^.Q is non- vanishing only if two gluons have 
positive and two of them have negative helicity. In eq. (|5.33|) , / and J denote the negative 
helicity gluons. 
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5.2 Results obtained using dimensional reduction 

The one-loop color subamplitudes obtained using dimensional reduction can most conve- 
niently be written with the help of the results given in the previous subsection. Obviously, 
to 0{e), only the divergent amplitudes may differ in the various schemes, therefore, we do 
not repeat the finite ones. 

The difference in the algebra when 't Hooft-Veltman or dimensional reduction schemes 
are used together with the tensor reduction of sect. 4, can be described by a single parameter 
which is chosen to be d-dimensional in the 't Hooft-Veltman scheme and has four- dimensional 
value in dimensional reduction. Changing the value of this parameter, one finds the following 
differences: 



xDf(-,±;T,+)-a^y(-,±;T,+) = cra4;o(-, ±; +) ( ^iV, - ^ ) . (5.34) 



c°f(-,±;T,+)-c«y(-,±;T,+) = crC4;o(-,±;T,+)^(iVc-^). (5.35) 

c?°f(-,±,T,+)-rf«y(-,±,T,+) = crrf4;o(-,±,T,+)^iVc, (5.36) 

while the m4;2, {m = a, c or d) subamplitudes are identical in the two schemes. We remark 
that our results for the gluon-gluon scattering without the fermionic contribution, obtained 
using dimensional reduction is identical with the corresponding expression obtained using 
string-based rules and the four- dimensional helicity scheme in ref. 

These differences may have two origins. They can arise from the different regularization 
of both the ultraviolet and the infrared singularities. The difference in the ultraviolet regu- 
larization can be absorbed into the change of Aqcd- This difference can be understood as a 
finite renormalization of the gauge coupling and can be calculated separately. In the back- 
ground field technique we obtain the finite renormalization from calculating the gluon self 
energy in two different regularization schemes and use the background field Ward identity 

= (Z:^)-^^ (5.37) 

The self energy has the well-known form. 



Klik) = ir (A;^A;, - k'g^,)U{k'), (5.38) 



where 



nHv,,^) . (|.)^ r(i . .) {-\ - 51^^) h- _ i , (5.39) 

in the 't Hooft-Veltman scheme and 

^ (_4^) V(l . e) (-1,3. - . ,5.40) 

in dimensional reduction. As a result of the finite renormalization, the difference between 
the two results vanish, therefore. 
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and 



(5.42) 



9 - Q 

After this renormalization, the difference between the two schemes for the amphtudes be- 
comes 

1 



fi^f(-,±;T,+) -a^|(-,±;T,+) = cra4;o(-,±;T,+) iVc 



c°f(-,±;T,+)-<r(-,±;T 



HV/ 



(5.43) 
(5.44) 



crd4-fi{-,±,T,+)^Nc. 



(5.45) 



The tilde in eqs. (|5. 43 - 5. 45 ) reminds us that the expansion parameter is the same in 
both schemes. In sect. 10, we shall discuss how these differences can be incorporated into 
the definition of the parton densities. 

As a final step in this section, following ref . , we calculate the change in Aqcd caused 
by the finite renormalization of the gauge coupling according to eq. (|5.42|) : 



A 



DR 



A^JgCXp 



The constant c is defined by 



1 +c- 



a 



MS 



2tt 



(5.46) 



(5.47) 



i.e., c = Nc/Q in agreement with the result of ref. obtained in Lorentz-Feynman gauge. 



6 Loop contribution to the next-to-leading order ma- 
trix element 

In this section, first we give explicit formulae how to obtain the loop contribution to the 
next-to-leading order matrix element summed over color and helicities and then compare 
the obtained results with the corresponding expressions given in ref. |Q. In this way, we can 
establish the transition rules for loop amplitudes which connect our results to those obtained 
in the conventional scheme. 

6.1 Squared matrix element in the 't Hooft-Veltman scheme 

It is useful to rewrite the F functions in terms of the physical variables s, t and m in a 
form which will be more convenient when the square of the matrix element is calculated. In 
the next-to-leading order amplitudes of sect. 5, we factored out the Born terms. Therefore, 
it is sufficient to present only the real part of the F functions since the imaginary part 
will not play any role in a next-to-leading order calculation of the squared matrix elements 
{TZeF = JF). To keep the crossing symmetry of the matrix element manifest, we leave the 
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sign of s, t and u open, and introduce an auxiliary mass variable Q to express the logarithms 

As mentioned previously, the analytic continuation is done according to formula 
To keep the real parts of the logarithms, we use the following substitutions: 



log 



Si4 

Sl2 
Sl3 



log(^) 



where 



£{x) = log 
One finds the following results. 
Process A: 



X 



log' (—) -^i2{s)+i2it)~2£is)m, 



log 



Si4 

2 f Sl2 
Sl3 



i2{s)+i2{u)-2i{s)i{u), 



log' r^) ^i,(u)+i2{t)~2i{u)m, 



and i2ix) = f{x) -7r'e{x). 



(6.1) 

(6.2) 
(6.3) 

(6.4) 



For the square of the matrix element only the sum of the Ta-i and J-'a-2 functions is 
needed. 



-^a s, t, u) = ne[F^^, {e, s, t, u) + F ^ {e, s, t, u)] 



2 \ £ 



V 



9 3\ in V /9 \ 9 



+ N, \-{l{u) - i{t)) + ^ + vr^ + 24(t) - 2e{t)e{u) - ^£{t) 

+ ^ (l - {m + 4(«) + vr^ - 2£(t)£(n)) + - £{u)) 

+ ^g(£(.)-£(n)) + 8 + 2^(£(u)-£(t)) 

s V s 

--/5o. 

e 



(6.5) 



(l2{t) + £2(n) + TT^ - 2£(t)£(?i)) - 4£(s)£(t) + Al{t)l{u) 



^{^^ s, t, u) = 7^e[Ft (e, s, t, u) + F 2^(£:, s, t, u) 



2 \ 



V 



\ e 
r2 



9 3\ in 1/ /9 \ 9 



1 



- m) + H + j2 + 2C2(() - 2<(()<(«) - ^e{t) 

-(£(s) - i{u)) + 8 - 4£(s)£(t) + A^{t)^{u) 



+- ( 1 - - ) U2{t) + £2(5) + TT^ - 2£(t)£(s)) + 2-(£(t) - £(s)) 



t 



e 



(6.6) 



15 



Process B: 

For the square of the matrix element, we shall need !F~^ in addition to the functions. 
T:^{e,s,t,u)^ (6.7) 



Process C: 



V_ 

1 



—.-A- -^Nf + - (-m - m + m] + TjA^/w 



9 



+ 



Nr. 



'-{l{s) - l{u)) + 8 - 2£(s)£(i) + 2£(i)£(M) 



Wr 



2e ~ 2 



+ 



2NJ{s)£{u) + 



(6.8) 



1 + - (i2{s) + i2{u) - 2e{s)i{u) + TT^) 



e 



\Nc\ 2e 2) "e" 



+ \ (J^(^) + \^^) - ^^2(5)) + iVc {^-Ji{u) + 4(5) + TT^ 



+ 



2Ar, 



^2(5) + ^2{u) + tt" + 2£(s)£(^/) + U{s) - U{u) 
+ ^ (1 + (1 - 2^) (£(.) - £(«)) 

+ (^)' (4(5) + 4(W) - 2£(5)£(W) + TT^) 

l + ^(l-(l-20 (£(.)- ^(x.)) 



e 



(6.9) 
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2 

e 

4 
+ - 

s 

1 

+ 7 

1 



^ (£(«) - m + ^ - ^{s)) + ^ - m 



It turns out that only the sum and difference of the J-'c-2 functions are needed. 
T+ie, s, t, u) = ne{F-+{e, s, t, u) + F+i{e, s, t,u))^ 

Wj 

i{s)i{t)-e{s)e{u) 

I (Ms) + £2{u) + TT^) - 3e{s)e{u) + 2e{s)e{t) - 2e{t)e{u) 

I (4(5) + i2{t) + n^) - 3i{s)e{t) + 2i{s)i{u) - 2i{t)i{u 



u 



2 

e 

1 

+ 7 



+- 



V Uj 



Process D: 



3£ 9 3 ^ ' 



/i 



2 \ e 



11 



67 



11 



+ Nf 



St 










1 - 


(- 


--) 



-(^--2j (^^{s) + ^2{t) + ^-2^{s) 
+Po- m - e{s)) - |iV;4 + ^2 W + 2e{s)e{t)) ] 



e 
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The square of the matrix element summed over helicities can be expressed in terms of 
the T functions. 

Process A: 



hel col \ J 



(6.14) 



Process B: 



hel col 

/ ^2 ^2 ^2 ^2 



Process C: 



(6.15) 



EE[-^*-^]NLO = 2/cr 



hel col 



v_ 



xi2V- 



^3 ^3 

(•^ctr(«> + ^^r(s, u, i)) + —J^'tis, t, u) + —J^'tis, u, t) 



(6.16) 



AT, 



(^-i+(s, t, «) + ^+-(s, t)) + (^+-(s, t, «) + ^-+(s, u, t) 
s^T~{s^ t, u) + t^J^^{s, t, u) — u'^J^^{s, t,u) I . 



Process D: 



E E[-^*-^]nlo = 2/cr2Ar2^ (^,(s, t, u) + Ta{t, u, s) + ^^(w, t)) , (6.17) 

hel col 



where 



§2 g4 

:rrf(s, u) = 2— J'^Ti"(s, w) + 2— J^rf-f (s, li, i) + ^ [^d^{u, t, s) + Tj^{t, u, s)) . (6.18) 
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6.2 Comparison with existing results 

The matrix elements given in the previous subsection differ in finite terms from those calcu- 
lated in the conventional regularization scheme by Ellis and Sexton The difference can 
easily be described if one rewrites the Ellis-Sexton matrix elements in the form given in ref. 

i- 

Following ref. we denote a certain cross section (e.g., inclusive one-jet cross section) 
by X. At next-to-leading order, X is a sum of two terms, 

X = X[2 ^ 2] +X[2 ^ 3], (6.19) 

where X[2 —>■ n] is the [2 n] part of the cross section. According to the factorization 
theorem, the physical cross section in the QCD improved parton model for hadron-hadron 
scattering is a folding between the parton densities and the hard-scattering cross section: 

X[2 n] = V / dXa f dXbfa/AiXa,lJ,)fb/B{Xb,(J')daa,b{XaPA,XbPB,fl,as{fl)). (6.20) 

In eq. ( |6.20|) , d&a^b is the hard-scattering cross section for the process a + & ^ ji + . . . + j„. 
It is defined as a product of the flux factor and the integral of the squared matrix element 
over the phase space of the final state particles 

daa,, = -, E (6-21) 
n! 2xaXbS 

where Sn{Pj-) is the so called measurement function that defines the infrared-safe physical 
quantity. The counting factors 1/n! are present when all partons are treated indistinguishable 
and we sum over the possible parton types. For the [2 — > 2] process, the square of the matrix 
element — summed over final spins and colors and averaged over initial spins and colors — 
has the following perturbative expansion 

(^IMia + b^j,+ = ^^^^^^^^ ^i^^'\d.p) + 2/ (^^^ cri^^'\a,p) + 0{g')Y (6.22) 

where we denote a = (a, 6,ji,j2), P = {PayPbyPji^P^) ^i^) represents the number of 
spin and color states of a parton type a. Comparing our notation to the one just defined, 
we see that 



2/ ^ cr^^'\3,p} ^J:J:[M*MUo. (6.23) 

\^ / hel col 

In ref. 0, using the results of Ellis and Sexton, the following structure has been found 
for the next-to-leading order term: 



^(6)(a,p) = 7/>(4)(a, 



P)|-iECK)-^E7(an)| (6.24) 



+ W^ogr-P^]^l^-\a,p) 

^ m<n \ V / 
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where ip^^\a,p) is the d-dimensional Born term, %jj^^\d,p) are the color-hnked Born squared 
matrix elements in d dimensions as defined in Appendix A of ref . and tjj^l (a, p) represents 
the remaining finite terms. The sum over m and n runs from one to four. In eq. ( |6.24|) , 



C{a) is the color charge of parton a and the constant 7(a) represents the contribution from 
virtual diagrams to the Altarelli-Parisi kernel. Specifically, 

C{g) = N,, 7(^?) = IPo, (6.25) 

V 3V 
Ciq) = 7(g) = (6.26) 

2iv,' '^^^ ^ ' 

Using this notation, the transition from the conventional scheme to the 't Hooft-Veltman 
scheme can be achieved by substituting the four- dimensional expressions for the ip'^^'' and 
'ip':^n^ functions and leaving the V'ns functions unchanged. In sect. 8, we shall analyze the 
structure of the squared matrix element for the Bremsstrahlung part of the cross section and 
shall point out that part of the cross section undergoes analogous change when passing from 
the conventional to the 't Hooft-Veltman scheme. 

Next, we consider the squared matrix element in the 't Hooft-Veltman and dimensional 
reduction schemes in the form as given in eq. ( |6.24| ). The difference described at amplitude 



level by eqs. ( ^.43| - p^ ) means a difference in the V^^s functions. We see, however, that the 
difference is proportional to the Born term. Conseqently, we can use the same (properly 
modified) ■j/'^fg functions in both schemes and include the difference into a scheme dependent 
modification of the 7(a) functions. Instead of 7(a), we use the following e-dependent 7(0,5) 
functions: 

7(a,£) = 7(a) +£7(a), (6.27) 

where 

l{g) = ^iV„ 7(g) = ic^, (6.28) 
6 2 

with Cp = V/{2Nc). We obtained the above expressions for the 7(a) terms from the differ- 
ences among the squared matrix elements in the 't Hooft-Veltman and dimensional reduction 
schemes for all processes. These differences determine an overconstrained system of linear 
equations (four equations for two unknowns). We could solve this overconstrained system 
consistently, which indicates that using the 7(0,5) functions in eq. ( p.24| ), universality will 
be maintained. We receive further indication of universality if we recall our observation that 
for the four-gluon amplitudes in pure gauge theory, one obtains identical results using the 
string-based rules and four- dimensional helicity scheme as in field theory using dimensional 
reduction. Extrapolating this rule to the five-gluon amphtudes, one can easily check, using 
the results of ref. 0, that universality is maintained at the five-point level. We have an 
additional confirmation that the 7(0,5) functions are process independent. We can obtain 
them if we extend the validity of the momentum sum-rule to the 5-dependent part of the 
Altarelli-Parisi kernels in the 't Hooft-Veltman scheme, as can easily be checked using the 
expressions of Appendix B. 

When we write the results for the squared matrix element in the form of eq. ( |6.24|) , 



but in terms of 7(0,5) functions, then the transition from the 't Hooft-Veltman scheme 
to dimensional reduction is accomplished by setting 7(a) = 0. This observation will be 
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important when we discuss how to shift the difference in the loop amphtudes described by 
7(a) into the parton density functions. 

7 Scattering amplitudes in N=l supersymmetric 
Yang-Mills theory at one loop 

The classical Lagrangian of the N=l supersymmetric Yang- Mills theory in Wess-Zumino 



gauge for the component fields reads pS 



C- = -\{F;.f-\>^rK (7.1) 

where F^^ is the usual Yang-Mills field strength of a vector field g in the adjoint repre- 
sentation, V is the usual covariant derivative and A is a Weyl spinor in the adjoint repre- 
sentation. From this Lagrangian one immediately sees that the kinematic structure of the 
gluon-gluino-gluino coupling in N=l supersymmetric Yang- Mills theory is the same as the 
gluon-quark-antiquark coupling in QCD. The only difference between these two couplings 
is the color charge: in QCD the quarks are in the fundamental representation while in the 
supersymmetric theory the gluinos are in the adjoint representation. Therefore, changing 
the color matrices in our previous calculation for the QCD amplitudes to F", we can 
obtain the gluino-gluino, gluon-gluino and gluon-gluon scattering amplitudes in N=l su- 
persymmetric Yang-Mills theory. Such a procedure has a direct application in tree-level 
QCD calculation: the subamplitudes for gluon-gluon scattering can be obtained from those 
of quark-gluon scattering ||2^, At one- loop level the color structure of the amplitudes in 
QCD does not allow for such a direct use of supersjnnmetry. One obvious application is going 
into the supersymmetry limit and use supersymmetry Ward identities as important checks 
on the calculation (see below). A more direct application is the following. Calculate the 
quark-gluon scattering amplitudes, change the color charges to obtain the gluino-gluon scat- 
tering amplitudes. Using supersymmetry Ward identities, obtain the gluon-gluon scattering 
amplitudes in N=l supersymmetric Yang-Mills theory. Calculate that part of the gluon- 
gluon scattering amplitude in QCD which is proportional to the number of quark flavors, Nj. 
This calculation is much simpler than the complete calculation because it involves the eval- 
uation of only those Feynman diagrams which contain a closed fermion loop. Subtract this 
contribution from the result in the supersymmetric theory with color charges corresponding 
to the adjoint representation of the fermions (i.e, Nf and simultaneously add it with 

the QCD color charges. Thus we obtain the correct amplitudes for gluon-gluon scattering in 
QCD saving considerable amount of work in that part of the calculation which is the most 
difficult to carry out. 

In order to obtain a meaningful result in a supersymmetric theory beyond tree level, one 
has to use a supersymmetric regulator or alternatively one has to restore the supersymmetry 
Ward identities if a supersymmetry breaking regularization has been used. At one-loop level, 
dimensional reduction is known to respect supersymmetry. It is, however, interesting to see, 
how supersymmetry can be restored when a non-supersymmetry-preserving regulator (such 
as the 't Hooft-Veltman scheme) is used. Therefore, we shall elaborate on both approaches. 
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Let us now give the results of the calculation performed in the dimensional reduction 
scheme. We shall use a self-evident extension of naming the processes, i.e., process B will 
be the gluino-gluino scattering, process C will be the gluon-gluino scattering and process 
D will be the gluon-gluon scattering. The results will be given for one gluino flavor. In a 
super symmetric theory, the amplitudes for those helicity configurations which are vanishing 
at tree level have to vanish at one loop as well. Explicit calculation shows that this property 
is fulfilled in the one-loop calculation. In the following, we give only those non-vanishing 
amplitudes from which one can obtain other amplitudes using parity transformation or cyclic 
property of the amplitudes. As usual, all particles are outgoing. The results can be given in 
a concise form: 

^USY,DR^_^ ±, +) = crmSusY(_^ ±^ +)N,F,-^{e, s^, ^is, ^m), (T.2) 
where m stands for b, c or d. The Fi'^ functions are universal: 

Ff\e,si2,si3,si4) = (7.3) 

4 3 2 si, o\ 3 
log 6 + n 

£ £ Si2 / e 

+ 6+' fs^logfii-^'^^'^^ 





\ Si3 Si2 2 Si3 

It is easy to check that these results are indeed supersymmetric, i.e., they satisfy certain 
on-shell supersymmetry Ward identities. These Ward identities can easily be derived noting 
that the supersymmetry charge Q{r]) with rj being the fermionic parameter of the trans- 
formation annihilates the vacuum. Then it follows that the commutator of Q{t]) with any 
string of operators creating or annihilating of a gluon or a gluino has a vanishing vacuum 



expectation value pBI. This statement is true in any order of perturbation theory. If Oj 



represent these operators, then we find the following supersymmetry Ward identities: 

= 



i=l . 



) =J2{ai...[Q,ai]...an)o. (7.4) 

«=i 

In the N=l supersymmetric Yang-Mills theory we considered above, the stand for 
and A^. The superscripts ± refer to the two possible helicity states of the vector and spinor. 



The action of the supersymmetry charge on the doublet {g, A) is as follows [pq] : 

[Q{^),gHp)] = TT^{p,v)^^- (7.5) 

[Q{rj),\^{p)] = Tr^{p,r^)g^. (7.6) 

Substituting these commutation relations into eq. ( |7.4| ) we obtain a relation among vari- 
ous scattering amplitudes for particles with different spin. These relations are the on-shell 
supersymmetry Ward identities referred to above. 

We remark that the supersymmetry Ward identities hold separately for each of the sub- 
amplitudes in which one can expand the full amplitude 



In order to find explicit relations, the V^^p,!]) functions are to be specified yet. V^^p,!]) 
is a complex function linear in the anticommuting c-number components of t] and satisfies 

r+(p,r7) = {T'{p,r])y = f]u-{p), (7.7) 
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where U-{p) is a negative helicity spinor satisfying the massless Dirac equation with momen- 
tum p. Using the freedom in choosing the supersymmetry parameter t], we choose it to be 
a negative hehcity spinor obeying the massless Dirac equation with arbitrary momentum k 
times a Grassmann variable 6. Then, 

r+(p, k) = r+(p, r]{k)) = 9{k + \p-). (7.8) 

As an application, we show how to derive the &f.o^^ amplitudes from the well-known 
results for ^4-0 (see eq. ( ^.33|) ). To obtain (Aj~A2 A3"A4 )o, consider {[Q,gig2g3^t])o 

= {[Q,9i 92 9t9t])o = T-ip,,k){X^g^gtXt)o (7.9) 
+r-(p2, k){g^X2gtK)o - r~(p4, k){g^g2gtgt)o, 



= ([g, Ar^?2"A^A+])o = -r-(p2,0(ArA2-A^A+)o (7.10) 

+r~(P3,0(Ar^^^3'A|)o - r"(p4,/)(Ar^2"A^^^)o. 

Choose k = p2 and I = p^. Using the properties of the spinor products and momentum 
conservation, we find 

(ArA2-A+A|)o = 64T^(-, -, +, +) = -iMM. (7.11) 

In an exactly analogous fashion, from {[Q , gi gt g^ "^t]) ^'^'^ ([Qj A^fi'^ A3 A4 ])o with k = p^ 
and I = p4, we obtain 

(ArA^ AjADo = blf-{-, +,-,+) = iMM^H. (7.12) 

■514 512 

We now see that our results given in eq. ( |7.2| ) are indeed supersymmetric due to the 
universality of the F^^ functions. When applying a supersymmetry transformation, the 
only change one finds is such that the Born terms transform into one another. 

From the results of sect. 5, we see that the color subamplitudes calculated in the 't Hooft- 
Veltman and dimensional reduction schemes only slightly differ. Consequently, one expects 
that the 't Hooft-Veltman scheme breaks supersymmetry only by a small amount at one-loop 
level. This is indeed what one finds in a direct calculation: 



m. 



SUSY.HV^ ^ j_ , ^ _ ,^SUSY 
4;1 



6„ 



-, T, ±, +) = crm|"^o' T, ±, +)N, Ff ^(e, Su, Sn, Su) - ^ • (7.13) 



with the same process independent Fi ^ functions as in eq. ( |7.2| ), but with process dependent 
5m- Sb = 5, 5c = 3, 5d = 1. 

According to the supersymmetry Ward identities (ff^), have we used a supersymmetry- 



preserving regulator, the constants 6m would have been process independent in eq. ( |7.13|) . 
In the followings, we shall argue that the supersymmetry breaking in the 't Hooft-Veltman 
scheme can consistently be restored applying a special form of finite renormalization. We 
obtain this renormalization via calculating the gluon and gluino self energies to 0{e) order 
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in 't Hooft-Veltman scheme and dimensional reduction scheme and require that the finite 
difference between the resuhs can be compensated by a finite renormahzation. For the self 
energies, we have the well-known expressions 

UfM = ^5''\k,K - fcV)n(A:'), (7.14) 
iE}(fc) = 5iE(P)^. (7.15) 

The functions n(fc^) and S(A;^) are regularization-scheme dependent. Explicit calculation 
for space-like k in the background-field gauge gives 



n"^M = [-^j T(l+e)(^----j+Zi-l + 0{e). (7.16) 

EHV(A,2) ^ ;v,f^)Y-i^)'r(l+E)fi + l)+Z{-l + 0(e), (7.18) 



e 



E°^P) = ^c(l^)' (-^y r(l + .)(! + 2) +0(.). (7.19) 

Requiring that the results are equal in the two schemes, we obtain the desired finite renor- 
malizations: 



and 

Zi = l + N,(^^y (7.21) 

So far, this renormahzation is very similar to that of sect. 6, where we discussed the dif- 
ference between Aqcd in the 't Hooft-Veltman and dimensional reduction schemes. There 
is, however, an important difference. The renormalizations discussed in the present section 
have to vanish as A;^ — *• 0. The reason is the following: In the spirit of dimensional regu- 
larization, the ultraviolet and infrared divergencies cancel exactly for on-shell propagation 
independently whether 't Hooft-Veltman scheme or dimensional reduction has been used. 
Therefore, the only breaking of supersjTiimetry could occur when calculating the loop cor- 
rection for off-shell lines. This implies that the renormahzation insertions occur on off-shell 
propagator lines only. Of course, in order not to break gauge invariance, additional finite 
renormahzation of the vertices is necessary according to the Ward identities 

Z{aa = Zi{Z{f/\ Z{,, = Zf{Z{y/'Zi (7.22) 

where Z^ is the gauge-coupling renormahzation^ From these equations we obtain the fol- 
lowing condition among the renormahzation constants 

Z{AAiZ{)-' = Z{,,iZl)-\ (7.23) 



^We note that in the helicity formalism the renormahzation of the four-gluon vertex at the one-loop level 
is irrelevant. 
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This condition is sufficient to ensure that if the renormahzation procedure is consistent for 
the four-point amphtudes then it is consistent for any n-point one-loop amphtude. The 
proof goes by induction. We assume that the finite renormahzation consistently restores 
supersymmetry for the {n — l)-point amplitudes. Attaching an external gluon leg on an 
(n — l)-point tree amplitude on a gluon line (either propagator or external leg) we introduce a 
finite renormahzation of Z^j^j^{Z^)~^ , while attaching the external gluon on a fermion line we 
introduce a renormahzation of Z^^^^{Zl)~^. According to eq. (|7.13|) , these renormalizations 
are equal, therefore, if the renormahzation was consistent at (n — l)-point level, it remains 
consistent at n-point level. Using this argument, and the results of ref. for the five- 
gluon amplitudes in the pure gauge sector presented in the 't Hooft-Veltman scheme, we can 
deduce the field theory result in dimensional reduction which turns out to be identical to the 
result in the four-dimensional helicity scheme. This is a further indication that dimensional 
reduction in field theory gives the same result as the four-dimensional helicity scheme when 
the string based rules are used. 

The consistency of our renormahzation procedure at the four-point level can be seen by 
performing the finite renormahzation on formula (|7.13| ) explicitly. We find that the constants 
6m in eq. ( |7.13| ) become universal, hence the required supersymmetry Ward identities are 
fulffiled. However, the absolute normalization is not set by condition (|7.23| ), consequently, 
the actual value of Sm is unknown. In order to set the absolute normalization, we require 
the background-gauge Ward identity between the gauge-field and gauge-coupling renormal- 
izations]^ ||2y], namely 

Z/ = (zi)-'/= = l-iiV.(|^)'. (7.24) 

Having the absolute normalization set, we obtain Sm = 0, i.e., we recover the result obtained 
using dimensional reduction exactly together with the correct shift in Aqcd between the two 
schemes. 



8 Bremsstrahlung contributions in 't Hooft-Veltman 
and dimensional reduction schemes 

In the ffist part of this paper, we have seen that one-loop amplitudes of massless QCD can 
be calculated using the 't Hooft-Veltman version of dimensional regularization and a simple 
rule for the transition to the conventional scheme has been found. The advantage of using 
the 't Hooft-Veltman scheme is obvious: both the calculation and the final result is much 
simpler than in the conventional scheme. We have also seen the usefulness of performing 
the loop calculation using a supersymmetry preserving or weakly violating regulator. On 
the other hand, the most widely accepted regularization scheme for the regularization of 
infrared divergencies which is also used in the factorization procedure is the conventional 
scheme. These facts motivate a thorough study of the various dimensional regularization 
schemes to establish the transition rules among them. We already know, how to make the 

^We use the background field formalism to the gluon sector only, therefore, we have to use the non- 
supersymmetric condition. 
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transition for the loop corrections. In order to establish the transition for physical quantities, 
we have to investigate the Bremsstrahlung contributions as well. 

To start with, we recall the form of the Bremsstrahlung contribution found in ref. In 
order to keep the discussion reasonably concise, we give the integrals in a schematic form 
and spell out those factors explicitly which depend on the regularization scheme. We admit 
that it is difficult to follow our discussion without consulting ref. for more details and 
precise definitions. However, we find it pointless to recapitulate every details which can be 
found in ref. 

The structure of the real corrections after integrating over the phase space is same as eq. 
( |6.24| ) for the virtual corrections, although the actual expression is far more complicated. 
Schematically we write 

4 

J[2 ^ 3] = ^ J[2 ^ 3]i, J[2 ^ 3]i = J[soft], + J[coll]i + J[NS]„ (8.1) 

i=l 

where X[2 3] is the [2 — ^ 3] cross section. To obtain X[2 — 3], we have to integrate over 
the momenta of the three final state partons. We can always arrange the calculation such 
that we label by 3 the parton with the smallest transverse momentum. The squared matrix 
element can be decomposed into a sum of four terms, where one term (|^5''®'^p)i {i = 1, 2, 3, 4) 
has singularities only when parton 3 is soft or coUinear to parton i. Then X[2 — > 3]j is the 
part of the cross section obtained by integrating (l^^^'^'^p),; X[soft]j, X[coll]i and X[NS]j are 
its soft, coUinear and finite terms respectively. The actual form of the X[soft]j, X[coll]i and 
X[NS]j integrals depends whether parton i is in the initial or final state. 

8.1 Soft contributions 

In the soft limit Ps — > 0, the variables of parton 3 can be integrated analytically. One is left 
with an integral over the variables of partons 1 and 2: 

X[soft].= J DMa,b)S2{pi,pl)ijf\3,p}, (8.2) 

<i,b,ji,j2 

where the sum runs over all possible fiavors of the initial (a and b) and final (ji and j2) state 
partons, Di is a ubiquitous prefactor which does not depend on the regularization scheme. 
The function L{a, h) describes the parton luminosity, S2{p^^,p^^) is the measurement function 
and the function ip^°^^{a,p) is the soft limit of the squared matrix element integrated over 
the variables of parton 3. They all depend on the regularization scheme. However, the 
dependence in the parton luminosity and the measurement function is irrelevant as far as 
our considerations are concerned because the same type of integral occurs when the [2 2] 
matrix elements are integrated over the phase space (see eqs. (29) and (33) in ref. P). 
When i is a final state parton, the function ■ipf^^{a^p) has the form 

^r'n«,P) = ^^'^{a,p)[^^C{a,)--2C{a,)\og{x)^ (8.3) 
+ E {-^log (^^) +i^m{x) + 0(.)| , {^ = 1,2), 
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while for the case when i is an inital state parton, we have 



^Pf\a,p) = ^W(g,p-)|lc(a.) + i2C(a,)log(^)} (8.4) 
+ ^ i^t\S,P) |-^log (^^) + 0{e)^ , {^ = 3,4). 

There is a term proportional to — {l/e) log(x) in eq. ( |8.3| ) (the same x appears in the argument 
of the lim functions too). In this term, x is an arbitrary positive number, less than one and 
it represents the fraction of pi that sets the upper limit of the integration in p^. It is used 
to separate the soft and collinear contributions in the soft-collinear region. The limit x ^ 1 
means that in the collinear region all momentum configuration is considered collinear except 
when p3 ~ 0. In ref. 0, x = 1/2 has been chosen for the sake of convenience. For our purpose 
of establishing the transition rules among the regularization schemes, it is more suitable to 
keep X arbitrary. The reason will be clear when the collinear integrals are investigated. 

Using the derivations in Appendix A of ref. |Q, we find that in eqs. ( ^.3] ) and ( ^.4] ) 
only the ip'^'^^ and iplm'^^ functions depend on the regularization scheme: for the conventional 
scheme they are the d- dimensional expressions, while for the 't Hooft-Veltman or dimensional 
reduction schemes they are to be taken in four dimensions. 

8.2 Collinear contributions 

The integral for the final state collinear singularities has an analogous form to that of the 
soft terms in eq. (|8^): 

X[coll],= j DMa,h)S2{p^,.p%W\^.p)^ (^ = 1,2). (8.5) 

In eq. 

is the collinear limit of the squared matrix element integrated over the variables of parton 
3. The function Z[a,x) is an integral of a finite expression: 

Z{a,x) = I dz(l — z)~ 

J X 

The lower limit of the integration is the same x we have met in eq. ( p.3|) . Pi)/a{z,e) is the 
usual Altarelli-Parisi splitting function in (4 — 26) dimensions without 2; — > 1 regulation. 
We shall term them e-dependent Altarelli-Parisi kernels. In Appendix B, we collected these 
e-dependent Altarelli-Parisi kernels in all three regularization schemes. Using those expres- 
sions, one can easily verify that Z{a,x) does not depend on the scheme in the limit x — 1. 
We are allowed to take that limit, because the sum of the soft and collinear terms is inde- 
pendent of X (as can be checked using the precise formulae of ref. 0]). Therefore, we find 



-2e 



2C{a) 



1.7) 



27 



again that in the final state colhnear term, eq. ( |8.5| ), apart from the irrelevant dependence 
in L and S2, it is only the Born term, ip^'^\d,p) which is regularization-scheme dependent. 

The structure of X[coll]j for initial state parton is more complicated. It can be split into 
two terms. The first has the structure encountered at the final state terms, while the second 
contains the the initial state collinear pole that is involved in the factorization theorem|^: 



X[coll]i 
In eq. (|878|) 



3,4). 



V^™'' has the form 



7(ai) - 2C(ai)log 
7(ai) - 2C{ai) log 



Xi 



(8.9) 



Xi 



Again, the whole regularization-scheme dependence in is contained in the Born term. 
The other term, G™^^ has a more complicated structure. As always, we emphasize only that 
part which depends on the regularization scheme. 



Gf\d,p) = -G 



(2-2) 
CT,j 



^.10) 



-a. 



2 P2 



P2 



2s2 \2nfx^ 
-l-finite terms. 



In eq. ( |8.10| ), G^^f^ is the part of the MS factorization counter term that is associated with 
particle i, 



G 



(2-2) 
CT,i 



(8.11) 



a. 



2P2 P2 



-2e 



^p^'Hd,p)S2{p^„p^jY: 



2s2 \27lfl 



tuia] 



uj{ai) Jx, z 



1 dz 



L \ a,b. 



Xi 



(4vr) 



z J eTil-e) 27r 



-p: 



ai/a\ 



It cancels exactly, when the factorization counter term is added. The Born term and the 
factors a; (a) which represent the number of spin and color states of a parton type a are 
four- dimensional expressions in the second line of eq. (|8.10|) . P' represents the e-dependent 
part of the £-dependent Altarelli-Parisi kernels. 



d ~ 



(8.12) 



The functions P' are given explicitly in Appendix B. The scheme dependence in eq. ( |8.10| ) 
occurs in the P' functions only. 

Summarizing this section, we conclude that the dependence on the regularization scheme 
in X[2 — *■ 3] appears only in the tp'^'^^ and ipl^'^^ functions and in a folding of the parton 

^°In ref. H, the whole contribution has been called 0^°'^ For our purposes this decomposition is useful. 
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luminosities with the e-dependent part of the e-dependent Altarelh-Parisi kernels, P' . The 
transition rule from the conventional scheme to the 't Hooft-Veltman or dimensional reduc- 
tion schemes is very simple: the ip'^^'> and i^l^^ functions has to be taken in four dimensions 
and L®P' has to be chosen in the proper scheme according to the formulas in Appendix B. 
This transition rule is in accordance with the transition rule found in sect. 6 for the virtual 
contributions. 



9 Transition rules for the hard-scattering cross section 
among dimensional regularization schemes 

According to the factorization theorem, the physical cross section in the QCD improved 
parton model for hadron-hadron scattering is a folding between the parton densities and the 
hard- scattering cross section: 

^(^{PA.Pb) =^ I dXa f dXbfa/A{Xa,fi)fb/B{Xb,fi)daa,b{XaPA,XbPB,fJ',as{fJ'))- (9.1) 
a,b 

In this equation, the fa/ a parton densities are process independent, their Q^-evolution is de- 
termined by the Altarelli-Parisi equations. At next-to-leading order the regularization and 
factorization scheme dependence of the Altarelli-Parisi kernels is cancelled by that of the 
hard-scattering cross section, daa,b- Therefore, in principle all schemes are equally accept- 
able. In practice, however, one is forced to use the existing parton density functions which are 
obtained by fitting a large set of data in a rather complicated phenomenological procedure. 
This fitting and the Q^-dependence of the fitted function is always worked out using the 
conventional regularization and MS (or DIS) factorization schemes. As a result, the existing 
parton density functions can be used only with those hard-scattering cross sections which 
are calculated in the conventional dimensional regularization scheme. On the other hand, 
as we have seen above, convenient applications of the helicity method to loop calculations 
require that regularization is performed in the 't Hooft-Veltman or dimensional reduction 
schemes. Therefore, it is vital to find the transition rules which tell us how to transform 
the expressions for the hard-scattering cross sections obtained in the 't Hooft-Veltman or di- 
mensional reduction schemes into the corresponding expressions in the conventional scheme. 
That goal is easily achievable using the results of sections 6 and 8. 

The cancellation of infrared divergencies can be seen manifestly when the analytic struc- 
ture of the infrared singularities are exhibited according to eqs. (|6.24| ) and ( |8.2| - |8.1i|) . In 



this cancellation the tp^^\d,p) and tp^'^\d,p) functions appear as formal objects, the only 
requirement is that they are the same in the real and virtual contributions and in the fac- 
torization counter term. We have seen that it is indeed the case in all three versions of 
dimensional regularization: in the conventional scheme, they are (i-dimensional expressions, 
while in the 't Hooft-Veltman and dimensional reduction schemes, they are taken in four 
dimensions. 

After cancellation of infrared divergencies and subtraction of initial-state collinear singu- 
larities, we find still process independent differences in the hard-scattering cross sections in 
the different regularization schemes. If we compare d&a^ in the 't Hooft-Veltman scheme to 
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that in the conventional scheme, then we find that only the P' functions are different. If we 
compare d(Ta,b in the dimensional reduction scheme to that in the conventional scheme, then 
we find that the difference is contained in part in the redefinition of Aqcd, in part in the P' 
functions and in part in 7(0, e): in the conventional scheme 7 takes its e-dependent value, 
while in dimensional reduction, its e-dependent part is set to zero. From these considerations 
we see that the only modification one has to perform on d(Ja,b obtained in non-conventional 
schemes in order to recover da^ {, in the conventional scheme is simply changing P' to the 
value in the conventional scheme, use e-dependent 7 parameters and the standard QCD 
parameter. 

In order to find the proper form of the loop contributions in the 't Hooft-Veltman or di- 
mensional reduction schemes which is useful for the cancellation mechanism (i.e., eq. (|6.24|) ), 



we have to calculate the ip^^\a,p) and i^mnK^iP) functions in four dimensions. But these are 
tree level calculations which can easily be done using the helicity method. For the definition 
of the 'ipmnK^yP) functions, we refer to ref. [Q. Knowing the ^/^(^^(a, p) and i'mnK^^P) func- 
tions in four dimensions and the complete loop contributions in d dimensions, we can deduce 
the form of the ^p^l{a,p) functions which constitute the finite part of the loop corrections 
in the conventional scheme. As far as the Bremsstrahlung contributions are concerned, in 
the 't Hooft-Veltman and dimensional reduction schemes, both the soft and the hard part 
of the cross section have to be evaluated in four dimensions. It is only the calculation of the 
P' functions in the 't Hooft-Veltman scheme which requires d-dimensional treatment (see 
Appendix B). 

We have given the transition terms and functions needed to recover the hard-scattering 
cross section in the conventional regularization scheme from calculations performed in the 
't Hooft-Veltman or dimensional reduction schemes. This is the main result of our paper, 
which was based upon the assumption that the analytic structure of the singularities found in 
ref. is process independent even after the minor modification obtained by the introduction 
of e-dependent 7 parameters. One expects that this structure is indeed universal because 
the singularities are linked to the external particles of the subprocess. To give more support 
to the universality, in sect. 10, we present further consistency check on the transition rules 
via going into the supersymmetric limit of QCD once again. 

10 Next-to-leading order Altarelli-Parisi kernels in the 
't Hooft-Veltman and dimensional reduction schemes 

We have seen that the only dependence on the regularization scheme in the hard-scattering 
cross section, daa,b — beyond the redefinition of Aqcd — is contained in P' and 7(a). We 
have established the transition rules to recover the hard-scattering cross section in the con- 
ventional MS scheme from a calculation performed in non-conventional schemes. Although, 
all phenomenological studies for the parton densities are performed using the conventional 
next-to-leading order Altarelli-Parisi kernels, it may still be useful to know these kernels 
also for dimensional reduction and 't Hooft-Veltman schemes. In particular, knowing the 
next-to-leading order Altarelli-Parisi kernels for a supersymmetry preserving regularization, 
such as dimensional reduction, has got its own interest. 
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The next-to-leading order Altarelh-Parisi kernels in dimensional reduction can be ob- 
tained in two different ways. On one hand, they can be determined in a direct calcula- 
tion using dimensional reduction. On the other, requiring that the physical cross section 
should not depend on the regularization scheme, the regularization-scheme dependence of 
the next-to-leading order Altarelli-Parisi kernels can be uniquely determined from that of 
the hard-scattering cross section. If the two approaches lead to the same result, then it 
would prove that unitarity is not violated as well as it would ensure the universality of our 
transition rules. A direct calculation of the next-to-leading order Altarelli-Parisi kernels in 
non-conventional schemes has not been carried out yet. Nevertheless, we can still perform 
for the dimensional reduction scheme a consistency test: we can test explicitly whether the 
next-to-leading order Altarelli-Parisi kernels satisfy the supersymmetry Ward identities for 
N=l supersymmetric Yang-Mills theory. 

For the sake of completeness, we review the derivation of the transformation rule for 
the next-to-leading order Altarelli-Parisi kernels using the example of deep-inelastic lepton- 
hadron scattering. The F2 structure function at next-to-leading order is 

-F2(x, Q) = E (u/a{x, /^) + ^ E fa/A{z, /i) ®. Cia{z, Q / ^^)\ + 0{al), (10.1) 
where the following notation has been introduced: 

f®,g=['dy dz f{z)g{y)6{x - yz) (10.2) 
Jo 

and i = q, q, a = q, g. In eq. QlO.lj ), Cia denotes the hard-scattering cross section — their 
actual form is irrelevant for our purposes]^. We can freely change Cia to C-^ = Cia + ^CioQ 
if we simultaneously change the parton densities to 

f'a/Ai^^f^) = fa/A{x,fl) - T^^Y. fb/A{f^) ®x ACabiQ/fi)l ■ (10.3) 



Indeed, writing eq. (|10.1|) in terms of C'ia and Z^/^, we find that the change in F2 is of the 
order 0{a^) which is neglected. The fi dependence of both fa/A and /^^^ is determined by 
the Altarelli-Parisi evolution equations, but with different kernels P and P' respectively: 



In the equation for /^/^, we substitute /^^^ according to eq. ( |10.3|) on the left hand side. 
Performing the differentiation and using 

= -ttPo— a., (10.5) 

we obtain 

E^a/.®./vA(/")= (10.6) 



E 



as ( A 



Pa/b ®x fb/A{f^) + ^ ( y/a/A ®x ACab " Y^{Pb/c ® Jc/a) ®x ^Cab 



+ 0{a 



In the language of deep-inelastic scattering, Cia is usually called coefRcient function. 
^^The DIS factorization scheme at next-to-leading order is defined by C^'^ = 0. 
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We can rewrite the right hand side in terms of /' and read off the change in the next-to- 
leading order AltareUi-Parisi kernel induced by a change in the hard-scattering cross section: 



a.. 



27r 



i/c 



0{al). (10.7) 



Table 1 contains the change in the hard-scattering cross section for making the transition 
from the conventional scheme to the 't Hooft-Veltman and dimensional reduction regular- 
ization schemes in the MS factorization scheme. These values can be obtained from the P' 
and 7(a) functions as described in Appendix B. Substituting these values for ACah into eq. 
( |10.7|) , we obtain the next-to- leading order Altarelli-Parisi kernels in the MS scheme. We 
list below only the difference in comparison to the standard MS expressions obtained in the 



conventional scheme 27 



Pfit^^) 



CFTjiNf{2 - -x-^ - -x^ + 2xlogx) + CgMI - x). (10.8) 



C^TfiiV/ (- 1 + a; + 2x log X -4xMoga;) (10.9) 

2 13 
CGTRNf{-l + -x~^ -4x+ —x^ - Axlogx + Ax'^ logx). 



pfii(^) 



p^M 



p^lU^) 



CpCG{2-'^x-'~^x^ + 2x\ogx) 



CFTRNf{-2 + \x-^ + \x^ - 2xlogx). 



10.10) 

;io.ii) 



PTM 



P^^i^) 



CpTRNf{3 - 



10 



-X 



X 



x^ + Ax log x) . 



= C pTriNf {-3 + 2x + x'^ -log X - Ax'^ log x) 

2 34 
+ CGTRNf{-2 + -x-' -10x + —x^ ~8x log x) 

+ T«iV;(7((7)-7(g))(l-2,T + 2x2). 



(10.12) 
(10.13) 



Cp(—x ""^ + -X — 2 logx — 2x logx) 



+ C^CG(2-a;~^ + x-2a;^ + 4a;logx) 

+ CF[-[3oX + {^{g)-^{q)){2-2x~'-x)]. 



(10.14) 



P^'i^) 



9 2 10 
= CFTRNf{-- + -x'^ -x + —x^ -Axlogx) 



(10.15) 



In these equations, we used the notation of ref. Cp = V/{2Nc), Cg = Nc and Tr = 1/2. 
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Now we turn to the N=l supersymmetric Yang-Mills theory discussed in sect. 7. In 
this case, if the regularization scheme is supersymmetric, the Altarelli-Parisi kernels are to 
satisfy the following identity to all orders in perturbation theory: 

P,i,i^) + P-oioi^) = P,d^) + Pmi^)- (10-16) 

The validity of this relation follows simply from the physical meaning of the Altarelli-Parisi 
kernels: they give the number density of the partons in the infinite momentum frame. There- 
fore, if we sum over the final states, it should not matter whether the initial state is a gluon 
or a gluino. The Altarelli-Parisi kernels are known to satisfy identity ( |10.16| ) to leading |27 
and in dimensional reduction to next-to-leading order 



Some care is to be taken when going to the supersymmetric limit of eqs. ( |10.12| - [T0.15|) . A 



virtual gluino loop on a gluon propagator contains a combinatorical factor of 1/2, while the 
gluon splitting into a gluino pair does not contain the same factor. Therefore, to obtain the 
supersymmetric limit, we take Tj^Nf —>■ in the splitting functions and 2TrA^j — > Nc in 
loop corrections. This rule can be written more consistently by using T^Nf N^-, (3q ^Nc, 
j{g) — > Nc/Q and 7(g) Nc/2 in addition to the usual rules, Cp —>■ Nc, Cq (this is 

the reason for not writing Pq, ^{g) and 7(g) in terms of Cg, Tr and Nf in eqs. ( |10.12| - [r0.15|) ). 
From eqs. ( |10.12| - [T0.15| ) for the combination 



A{z) ^ P,,,{z) + P~g,,{z) - Pg/^{z) - Pg/^{z), (10.17) 

we obtain 

ADR(^) _ aMS(^) = _ _ ^ 7^2 _ (2x + - 1) logx + -x-') (10.18) 

V6 3 3 / 

where Pa/h{z) is the Altarelli-Parisi kernel without the part that is proportional to 5(1 — 2;). 
Comparing eq. (|10.18|) to formula (19) of ref. |27|, we conclude that in the dimensional 



reduction scheme A (2;) = 0. This is the consistency check on our transition rules we re- 
ferred to above. This test, however, still does not prove that unitarity is maintained in the 
't Hooft-Veltman and dimensional reduction schemes. It may be violated by terms which 
are proportional to {N^ — Nf). 

Finally, we remark that our derivation of eq. (|10.18 ) is different from that of Antoniadis 



and Floratos ||2^ and can be considered as an independent confirmation of the result that 



in dimensional reduction, relation (|10.16|) remains valid in next-to-leading order. 



11 Summary 

We calculated the one loop corrections to the helicity amplitudes of all 2 2 parton scatter- 
ing processes in QCD using dimensional reduction as well as 't Hooft-Veltman regularization 
schemes for regularizing both the ultraviolet and infrared singularities. We explicitly demon- 
strated the cancellation of the soft and collinear singularities between the loop corrections 
and Bremsstrahlung contributions using the general structure of the Bremsstrahlung singu- 
larities described in ref. Using somewhat heuristic arguments concerning the differences 
between the results obtained in the 't Hooft-Veltman scheme and in the conventional di- 
mensional regularization scheme, we could reproduce the results of Ellis and Sexton for the 
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one loop contributions to the spin and color summed cross-sections of all 2 — 2 parton 
subprocesses. Using the general form of the soft and collinear singularities, we have found 
universal (process independent) transition rules which can be used to transform the results 
obtained in the 't Hooft-Veltman scheme and in dimensional reduction into those of the 
conventional regularization scheme. These transformation rules have practical importance 
because the phenomenological parton density functions determined from fits to the data 
assume next-to-leading order Q^-evolution as given in the conventional dimensional regular- 
ization scheme. Therefore, when a physical cross section is calculated, the hard scattering 
parton cross sections have to be given as obtained in the conventional scheme. Since we 
know how to translate the results between the above mentioned three schemes, one can do 
the actual calculation in the most convenient scheme. We note that in the 't Hooft-Veltman 
scheme and in dimensional reduction one can use the helicity method directly. 

With explicit calculation we demonstrated that in dimensional reduction, as expected, 
the supersymmetry Ward identities of the scattering amplitudes of the 2 — > 2 parton sub- 
processes of the N=l supersymmetric Yang-Mills theory are satisfied also in next-to-leading 
order. QCD differs from the N=l supersymmetric Yang-Mills theory only in the color 
representation of the quarks and gluinos, therefore, the supersymmetry Ward identities can 
be used to obtain a non-trivial test of the correctness of the calculation and to obtain the 
amplitudes of different subprocesses with the help of the supersymmetry Ward identities. 

We have shown that the supersymmetry Ward identities can be derived also using the non- 
super symmetric 't Hooft-Veltman regularization scheme by working out simple finite renor- 
malization factors. Using this result we find that the string theory based four- dimensional 
helicity scheme of Bern and Kosower gives the same result for the four- and five-gluon am- 
plitudes as the dimensional reduction scheme. 

It is not obvious that in the 't Hooft-Veltman scheme and in the dimensional reduction 
unitarity is not violated by some finite terms when the soft and collinear singularities are 
cancelled between the loop corrections and the Bremstrahlung contributions. This difficulty 
does not arise if we know the transition rules to the conventional regularization scheme. 
Using the transition rules we calculated the next-to-leading order Altarelli-Parisi kernels for 
dimensional reduction and for the 't Hooft-Veltman scheme. Unitarity is maintained in this 
case by construction. We have found that the next-to-leading order Altarelli-Parisi kernels 
of the N=l supersymmetric Yang-Mills theory obtained in dimensional reduction satisfies 
the supersymmetry Ward identities . 
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Appendix A 



In this Appendix, we spell out eq. ( [4 .21) explicitly for = 2, 3, 4 in a form which is used 
in the actual calculation. We also calculate all the necessary integrals. 
Case N = 2: 

-^ = 77^(^1 E(-p')'r(-A; + £) / dxidx2 5(l-xi-X2)(xiX2)'=-^Ar,(J,p), 

(A.l) 

where 

Ji = X2P, (A.2) 
J2 = ~xip. (A. 3) 

Nk{J,p) is a polynomial of Xi and X2- The general form of a term in the integral is 

/ dxidx25{l-Xi-X2)x'l^X2^{xiX2)''~'' = B{ni + k + l-e,n2 + k + l-e). (A.4) 

J [0,1] 2 

Case A^ = 3: 

D{x,p) = -{xiX2pI + X2XspI + x^XipI). (A. 5) 

For the 0—^4 partons case, two of the external legs are on shell. If i and j are the indices 
of the propagators joining to the massive external leg (pj), then eq. (^4.2| ) becomes 

,2\e 1 

(47r)2 \^ -p 

j ^ ^ ^ dxi dxj dxk 5{1 - Xi- Xj - Xk){xiXjY'^'^ Nk{J,p), 



M = -7^2 ] E(-p])'-'r(l -k + e) (A.6) 



where 



Ji = XjPj + Xk{pj+Pk), ( A . 7) 

Jj = XkPk + Xi{pk + Pi), (A. 8) 

Jk = XiPi + Xj{pi + Pj). (A. 9) 
Nk{J,p) is a polynomial of Xi, xj and Xk- The general form of a term in the integral is 

dxi dxj dxk (5(1 — Xi — Xj — Xk)x'^*x^j^ x^.'^ {xiXj)^~^~^ . (A. 10) 

[0,1] 3 

After changing variables. 



Xk = {l-x){l-y), 
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the integral becomes 

f^dx di/x"'(l-x)"^+"'=?/"^(l -?/)"'= (x(l - -x) = (A.ll) 

JO Jo 

B{ni + k — e^Uj + rik + k + 1 — e)B{nj + A; — e, + 1). (A. 12) 

Case iV = 4: 

D{x,p) = -{xiX2pI + X2X3PI + x^Xipl + X4X1PI + X1X3S12 + X2X4S14), (A.13) 

where 

Sij = {pi+Pjf. (A. 14) 
For the — > 4-parton case, the external legs are on shell. Therefore, eq. ( [4 .21) becomes 

^ = -(^2i'^''yT.n'2-k + e) (A.15) 



k=0 



/ dx,dx2dx,dx,6{l-x,-X2-x,-x,)- 

JfO,114 



Nk{J,p) 



/[0,1]4 {-XiXsSu - X2X4Sur-''+^ 

where 

^1 = X2P1 + X2,{pi + P2) + 3:4(^1 + P2 + Ps), (A.16) 

J2 = X3P2 + Xi{p2 +P3) + Xi{p2 +P3 +P4), (A.17) 

-''s = 2:4^3 + a;i(p3 +P4) +2:2(^3 +P4 + Pi), (A. 18) 

Ja = XiP4 + X2{Pa+Pi) + X3{Pa+Pi +P2)- (A.19) 

^k{JiP) is a polynomial of Xi, X2, X3 and X4. The generic form of the integral is 

l^[P{xi)] = (47r)-2+T(2 -k + e) (A.20) 



X^-/ dx,dX2dX3dX4^(l-X,-X2-X3-X4)- 



/[o,i]4 " " ' " ' " " " (,-a;ia;3Si2 - X2a;4Si4)2 

It is remarkable that the same type of integrals, but only for k = 0, appear in the 
calculation of one-loop amplitudes using the string-based approach [|l|. In ref. [^], a gen- 
eral method has been described how to calculate a large class of these one-loop Feynman 
parametric integrals. The main idea of the evaluation is derivation of a differential relation 
between the basic scalar integral (X'^[l]) and the so called scaled integrals which are 
related to X^[P{xi)] according to 

i'[Pixi)] = |^n«,j J^-[P(x,M)], (A.21) 

where are positive, real parameters, satisfying aia^ = (—512)"^ and a2a4 = (—514)"^. 

The basic scalar integral can be evaluated either using the partial differential equation 
technique of ref. |29|, or by direct integration. We have chosen the latter approach and 
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found agreement with the result of ref. [0]. For the cases k ^ the same method with trivial 
modifications works. 



Appendix B 



This Appendix is devoted to the derivation of the changes in the coefficient functions 
appearing in Table 1. 

The first step in the determination of the ACab functions is the derivation of the e- 
dependent Altarelli-Parisi kernels in all three regularization schemes. This accounts for 
finding the ^-dependent splitting functions without 2; — >■ 1 regulation. These functions in 
the dimensional reduction scheme coincide with the four-dimensional Altarelli-Parisi kernels 
first obtained by Altarelli and Parisi. We recall them to set the notation: 



P, 



9/g 



z 



2C, 



G 



Z 



1-z 



TrNj 



1 - z 

2 



zil 



z' + il-z) 



p, 



z 



1 + (1 



1 + z- 



(B.l) 
(B.2) 

(B.3) 
(B.4) 



We denote the e-dependent Altarelli-Parisi kernels by the regularization scheme-dependent 
functions Pa/b{z^ e) and decompose them into a scheme-independent and a scheme-dependent 
part: 

Pa,b{z,e) = Pa,b{z)+eP',„{z). (B.5) 

Our goal is to find the scheme-dependent functions, P'^i^{z). 

Clearly, when the initial parton is a quark, the only dependence on the dimension of space- 
time enters in the integral over the loop momentum. As a result, the splitting functions for 
quark splitting are the same in the 't Hooft-Veltman and conventional schemes. Repeating 
the calculation of Altarelli and Parisi in (4 — 2e) dimensions, one finds |31 



p',t\^) 



-Cpz, 
-Cf(1 - z) 



(B.6) 
(B.7) 



We follow the procedure of Appendix A of ref. to derive the splitting functions for 
gluon splitting. In that reference, the splitting for a polarized gluon were obtained in the 
form 



q/9^ 

9/9^ 



1 
2 

2a 



6'^ - 4z{l - z) 



G 



I- Z 



+ 



5'^ + 2(1 - e)z{l - z) 



(B.8) 
(B.9) 
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In eqs. ( pj.8D , ( |B.9| ), q is the transverse momentum of the product particle, and are the 
two components of q. The £-dependent AltareUi-Parisi kernel is then (see ref. 0) 



Pa/,{z,e) = ^^5,,Pl]^{z,c^). (B.IO) 



Both uj{g) and the contraction of the transverse indices are scheme dependent. In the 
't Hooft-Veltman scheme they are taken in four dimensions, while in conventional dimensional 
regularization they are in (i-dimensions. Using eqs. (|B.8[ - P.10[) , we obtain 



^r^W = 0, (B.ll^ 



KTi^) = -2T^Nfz{l-z) (B.12) 



in the conventional scheme, in agreement with the results of ref. [^, while in the 't Hooft- 
Veltman scheme, we find 



P^«7(^) = -2Cgz{1-z), (B.13) 
= 0, (B.14) 



9/9 
3/HV/ 
q/9 y 



in agreement with the results of ref. 0. 

The difference between the hard scattering cross sections when changing the regular- 
ization scheme determines the ACab{z) functions. These differences were analyzed in great 
details in the main text. The differences can entirely be described by the scheme dependence 
of the P' functions and by that of the 7(a) constants. The 7(a) constants represent the con- 
tribution from the virtual graphs to the AltareUi-Parisi kernel, therefore, we may interpret 
the 7(a) constants as the contribution from the virtual graphs to the e-dependent part of 
the ^-dependent AltareUi-Parisi kernel. The 7(a) constants are scheme dependent: 

^CDR(^)^^HV(^) = In,, (B.15) 

o 

^CDR(^) = ^HV(^) = Ic,, (B.16) 

7°^(^7)=0,7°%) = 0. (B.17) 
Using eqs. (^, ([B^), ( IBTIHB:!^ ), (|KT31 - PT7|) , it is now easy to fiU the entries of Table 1. 
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